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Charged particle assisted nuclear reactions in solid state environment: 
renaissance of low energy nuclear physics 
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The features of electron assisted neutron exchange processes in crystalline solids are surveyed. It 
is found that, contrary to expectations, the cross section of these processes may reach an observable 
magnitude even in the very low energy case because of the extremely huge increment caused by the 
Coulomb factor of the electron assisted processes and by the effect of the crystal-lattice. The features 
of electron assisted heavy charged particle exchange processes, electron assisted nuclear capture 
processes and heavy charged particle assisted nuclear processes are also overviewed. Experimental 
observations, which may be related to our theoretical findings, are dealt with. A possible explanation 
of observations by Fleischmann and Pons is presented. The possibility of the phenomenon of nuclear 
transmutation is qualitatively explained with the aid of usual and charged particle assisted reactions. 
The electron assisted neutron exchange processes in pure Ni, Pd and Li — Ni composite systems (in 
the Rossi-type E-Cat) are analyzed and it is concluded that the electron assisted neutron exchange 
reactions in pure Ni and Li — Ni composite systems may be responsible for recent experimental 
observations. 
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I. INTRODUCTION 

Since the ’’cold fusion” publication by Fleischmann 
and Pons in 1989 jlj a new field of experimental physics 
has emerged. Although even the possibility of the phe¬ 
nomenon of nuclear fusion at low energies is in doubt in 
mainstream physics, the quest for low-energy nuclear re¬ 
actions (LENR) flourished and hundreds of publications 
(mostly experimental) have been devoted to various as¬ 
pects of the problem. (For the summary of experimen¬ 
tal observations, the theoretical efforts, and background 
events see e.g. Q, 0.) The main reasons for revulsion 
against the topic according to standard nuclear physics 
have been: (a) due to the Coulomb repulsion no nuclear 
reaction should take place at energies corresponding to 
room temperature, (b) the observed extra heat attributed 
to nuclear reactions is not accompanied by the nuclear 
end products expected from hot fusion experiences, (c) 
traces of nuclear transmutations were also observed, that 
considering the repulsive Coulomb interaction is an even 
more inexplicable fact at these energies. 

Motivated by the observations in the above field we 
search for physical phenomena that may have modify¬ 
ing effect on nuclear reactions in solid state environment. 
Earlier we theoretically found a, a that if the reaction 
p+d —► 3 He takes place in solid material then the nuclear 
energy is mostly taken away by an electron of the environ¬ 
ment instead of the emission of a 7 photon, a result that 
calls the attention to the possible role of electrons. Con¬ 
cerning the assistance of the electrons and other charged 


constituents of the solid, a family of electron assisted 
nuclear reactions, especially the electron assisted neu¬ 
tron exchange process, furthermore the electron assisted 
nuclear capture process and the heavy charged particle 
assisted nuclear processes were discussed mostly in crys¬ 
talline solid state (particularly in metal) environment [6|, 
[7]. The aim of this paper is to summarize our theoretical 
findings and on this basis to explain some experimental 
observations. 

We adopt the approach standard in nuclear physics 
when describing the cross section of nuclear reactions. 
Accordingly, heavy, charged particles j and k of like pos¬ 
itive charge of charge numbers Zj and Zk need consid¬ 
erable amount of relative kinetic energy E determined 
by the height of the Coulomb barrier in order to let the 
probability of their nuclear interaction have significant 
value. The cross section of such a process can be derived 
applying the Coulomb solution <^(r), 

^(r)=e lk ' 1 7(k,r)/v / K, (1) 

which is the wave function of a free particle of charge 
number Zj in a repulsive Coulomb field of charge number 
Zk 0; i n the description of relative motion of projectile 
and target. In (JT]) V denotes the volume of normalization, 
r is the relative coordinate of the two particles, k is the 
wave number vector in their relative motion and 

/( k,r) = e —W 2 r(l + ipj k )iFi(-ir] jk , 1 ;i[kr - k • r]), 

( 2 ) 

where \F\ is the confluent hypergeometric function and T 
is the Gamma function. Since ip(r) ~ e~ 7rVjk//2 r(l+ir/j k ), 
the cross section of the process is proportional to 

2irr] jk ( E) 
exp [2nr]jk (E)] - 1 
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the so-called Coulomb factor. Here 


Vjk (E) = ZjZkOCf 


/ rn 0 c 2 

i aik ~2E 


(4) 


is the Sommerfeld parameter in the case of colliding parti¬ 
cles of mass numbers Aj , Ak and rest masses rrij = Ajuiq, 
rrik = AkiriQ. m^c 2 = 931.494 MeV is the atomic energy 
unit, a/ is the fine structure constant and E is taken in 
the center of mass (CM) coordinate system. 


Ujk — 


AjA k 

Aj + Ak 


(5) 


is the reduced mass number of particles j and k of mass 
numbers Aj and A Thus the rate of the nuclear reaction 
of heavy, charged particles of like positive charge becomes 
very small at low energies as a consequence of Fjk(E) 
being very small. 

In the processes investigated the Coulomb and the 
strong interactions play crucial role. The interaction 
Hamiltonian Hj comprises the Coulomb interaction po¬ 
tential Vcb with the charged constituents of surroundings 
(solid) and the interaction potential Vst of the strong in¬ 
teraction: 


Hi = Vcb + Vst ■ ( 6 ) 

(The Coulomb interaction between the charged partici¬ 
pants of the nuclear reaction is taken into account using 
©■) Therefore the charged particle assisted nuclear re¬ 
actions are at least second order in terms of standard 
perturbation calculation. According to (|6|, the lowest 
order of S-matrix element of a charged particle assisted 
nuclear reaction has two terms which can be visualized 
with the aid of two graphs. However, the contribution by 
the term, in which Vst according to chronological order 
precedes Vcb, is negligible because of the smallness of the 
Coulomb factor the root square of which is appearing in 
the matrix element of Vst in this case. (In the following 
we only depicts the graph of the dominant term.) 

When describing the effect of the Coulomb interaction 
between the nucleus of charge number Z and a slow elec¬ 
tron one can also use Coulomb function, consequently, 
the cross section of the process to be investigated is pro¬ 
portional to 


Fe(E) = 


2nr] e (E) 
exp [2nij e (E)} 


V 


but with 


Ve = —ZdLf 



(7) 

( 8 ) 


Here m e is the rest mass of the electron. In the case of 
low (less than 1 keV) kinetic energy of the electron F e (E) 
reads approximately as F e (E) = \2Tir] e (A 1 )! > 1. 

For instance, the cross section of electron assisted neu¬ 
tron exchange process (as it will be discussed later, and 


the graph of which is depicted in Fig. 1) is proportional 
to F e (E) only (instead of Fjk(E)) since the neutron takes 
part in strong interaction and so the corresponding ma¬ 
trix element does not contain Coulomb factor. The in¬ 
crement in the cross section due to changing Fjk(E) for 
F e (E ) in the case of electron assisted neutron exchange 
process can be characterized by the ratio F e (E)/Fjk(E) 
which is an extremely large number. The cross section 
of electron assisted neutron exchange process has a fur¬ 
ther (about a factor 10 22 ) increase due to the presence 
of the lattice since the cross section is also proportional 
to l/v c . Here v c ~ d 3 is the volume of the elementary 
cell of the solid with d the lattice parameter of order of 
magnitude of 10 -8 cm. The extremely huge increment in 
the Coulomb factor increased further by the effect of the 
lattice makes it possible that the cross section of electron 
assisted neutron exchange process may reach an observ¬ 
able magnitude even in the very low energy case. Thus 
it can be concluded that the actual Coulomb factors are 
the clue to the charged particles assisted nuclear reac¬ 
tions and therefore we focus our attention to them espe¬ 
cially concerning the Coulomb factors of heavy charged 
particles. 

It is worth mentioning, that usual nuclear experiments, 
in which nuclear reactions of heavy charged particles are 
investigated, are usually devised taking into account the 
hindering effect of Coulomb repulsion. Consequently, the 
beam energy is taken to be appropriately high to reach 
the energy domain where the cross section of the pro¬ 
cesses becomes appropriately large. Therefore in an or¬ 
dinary nuclear experiment the role of charged particle 
assisted reactions is not essential. 


II. APPLIED METHOD PRESENTED IN 
ELECTRON ASSISTED NEUTRON EXCHANGE 
PROCESS 

Recognizing the possibility and advantage of the assis¬ 
tance of electrons in LENR we consider first the electron 
assisted neutron exchange process, namely the 

e + %X + %Y -> e' + + £ +1 F + A (9) 

reaction jj] (see Fig.l). Here e and e' denote electron 
and A is the energy of the reaction, i.e. the difference 

between the rest energies of initial (jz]X +%* yJ and 

final (i 1 1 ~ 1 X+ 4 2+ly ) states. 

In © the electron (particle 1) Coulomb interacts with 
the nucleus X (particle 2). A scattered electron (par¬ 
ticle 1'), the intermediate neutron (particle 3) and the 
nucleus (particle 2') are created due to this inter¬ 

action. The intermediate neutron (particle 3) is captured 
due to the strong interaction by the nucleus ^ Y (particle 
4) forming the nucleus ^ +1 F (particle 5) in this man¬ 
ner. All told, in © the nucleus X (particle 2) looses a 
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FIG. 1: The graph of electron assisted neutron exchange pro¬ 
cess. Particle 1 (and 1’) is an electron, particle 2 is a nucleus 
which looses a neutron and becomes particle 2’. Particle 3 
is an intermediate neutron. Particle 4 is the nucleus which 
absorbs the neutron and becomes particle 5. The filled dot 
denotes Coulomb-interaction and the open circle denotes nu¬ 
clear (strong) interaction. 


neutron which is taken up by the nucleus z*Y (particle 
4). The process is energetically forbidden if A < 0. It 
was found, as it will be seen later, that the electron takes 
away negligible energy. In this process the Coulomb fac¬ 
tor of electrons arises only since the particle, which is 
exchanged, is a neutron. 

The physical background to the virtual neutron strip¬ 
ping due to the Coulomb interaction is worth mention¬ 
ing. The attractive Coulomb interaction acts between 
the Z\ protons and the electron. The neutrons do not 
feel Coulomb interaction. So one can say that in fact 
the nucleus ^~ x X is stripped of the neutron due to the 
Coulomb attraction. 

As an example we take Ni and Pd as target material. 
It is thought that the metal (Ni or Pd) is irradiated with 
slow, free electrons. In this case reaction © reads as 

e+ z ±x + z 2 X ->e'+ z^X-h % 2+1 X + A (10) 
with Z = Z\ = Z 2 . 

Now we demonstrate our calculation. Let us take a 
solid (in our case a metal) which is irradiated by a mo- 
noenergetic beam of slow, free electrons. The correspond¬ 
ing sub-system Hamiltonians are H so ud and H e . It is 
supposed that their eigenvalue problems are solved, and 
the complete set of the eigenvectors of the two indepen¬ 
dent systems are known. The interaction between them 
is the Coulomb interaction of potential V Cb (x) and the 
other interaction that is taken into account between the 
nucleons of the solid is the strong interaction potential 
V St (x) . In the second order process investigated an elec¬ 
tron takes part in a Coulomb scattering with an atomic 
nucleus of the solid. In the intermediate state a virtual 
free neutron n is created which is captured due to the 
strong interaction with some other nucleus of the solid. 
The reaction energy A is shared between the quasi-free 
final electron and the two final nuclei which take part in 
the process. Since the aim of this paper is to show the 
fundamentals of the main effect, the simplest description 
is chosen. 

The electron of charge —e and the nucleus X of 
charge Ze take part in Coulomb-interaction. We use a 


screened Coulomb potential of the form 

yCb M = -f^2 / ^72 ex P ( zf i • X ) rfc l ( n ) 

with screening parameter A and coupling strength e 2 = 
ctfhc. For the strong interaction the interaction potential 

r st (x) = ~/ exp t | ~( |x|) (12) 

is applied, where the strong coupling strength / = 0.08/Ic 
[9] and 1/s is the range of the strong interaction, (h is 
the reduced Planck constant, c is the velocity of light and 
e is the elementary charge.) 

According to the standard perturbation theory of 
quantum mechanics the transition probability per unit 
time (' Wfi ) of this second order process can be written 
as 

^ = ( 13 ) 


with 




v f^ b 

A 


(14) 


Hereof is the matrix element of the Coulomb potential 
between the initial and intermediate states and Vj* is the 
matrix element of the potential of the strong interaction 
between the intermediate and final states, furthermore 


AE^ = E tl -E i - A ifl . (15) 


Ei 1 E^ and Ef are the kinetic energies in the initial, in¬ 
termediate and final states, respectively, A is the reaction 
energy, and A ifJi is the difference between the rest energies 


of the initial 



and intermediate 


i^ 1 1 X and rt)j 


states. 


A = A_ + A+, A ifl = A_ - A n (16) 


with 


A_ = A Al - A j4i -i and A + = Aa 2 - A^+i- (17) 

A/tj, Aa 2 , A A2+ i and A„ are the energy excesses 

of the neutral atoms of mass numbers Ai, A\ — 1, A 2 , 
A 2 + 1 and the neutron, respectively. [10|. The sum 
of initial kinetic energies (Ei) is neglected in the energy 
Dirac-delta 5(Ef — Ei — A) and A E^ further on. 

Now for the sake of simplicity we reindex the particles. 
Particle indexed with e is the electron, particle indexed 
with 1 is initially the nucleus (J 1 X (particle 2 in Fig. 1) 
and finally ^ 1 ~ 1 X (particle 2’ in Fig. 1), particle indexed 
with 2 is initially the nucleus ^ 2 X (particle 4 in Fig. 1) 
and finally ^ 2+1 X (particle 5 in Fig. 1). 

E f = E fe (k /e ) + E fl (k,) + E f2 (k 2 ), (18) 
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E/j. — Ef e (k f e ) + E^i (ki) + E n (k„), (19) 


where 


E „ (K) = ^ ( 20 ) 

is the kinetic energy, k fj = k j is the wave vector and 
rrij is the rest mass of particle indexed with j in the final 
state (j = 1 , 2 ). 


E n (k„) 


2 ui n 


( 21 ) 


is the kinetic energy, k n is the wave vector in the inter¬ 
mediate state and m n is the rest mass of the neutron. 
Epi (ki) is the kinetic energy of the first particle in the 
intermediate state, and E^i (ki) = Ef i (ki). The kinetic 
energy of the electron in the initial and final state 


E iK = 


fr 2k ?e 

2m P 


and Ef e = 


h 2 k 2 


fe 


2 TOp 


( 22 ) 


with k i e and k/ e denoting the wave vector of the electron 
in the initial and final state. The initial wave vectors k^i 
and k ,2 of particles 1 and 2 are neglected. The initial, 
intermediate and final states are determined in Appendix 
A., the V°\ Vf* matrix-elements are calculated in Ap¬ 
pendix B. and the transition probability per unit time 
is calculated in Appendix C.. (Appendix D. is devoted 
to the approximations, identities and relations which are 
used in the calculation of the cross section.) 


A. Cross section of electron assisted neutron 
exchange process 

The cross section a of the process can be obtained from 
the transition probability per unit time (1751) dividing it 
by the flux v e /V of the incoming electron where v e is 
the velocity of the electron. 

c ajhcZ 2 E; 2 =-m 2 1*2 (k 2 )] 2 (23) 

Ve tt 3 v c (|kr + k 2 | 2 + A 2 )' (A£^)i n=kj 

( |Fl (k2)|2 ) A * rA > s ( E f ~ A ) rf3fc a, 

where v c is the volume of elementary cell in the solid, ta 2 
is the relative natural abundance of atoms A, 

(k 2 ) = J $ii (r„i) e~ tk2 A i- 1 r " 1 d 3 r ri i, (24) 

li=mi 

— £ |F l( k„)f (25) 

ll= — TOl 



and 


F 2 ( k 2 ) = J $} 2 (r„ 2 )e ik2 r " 


(26) 


-/ 


exp(-s^ i r„ 2 ' 


drr n2 . 


Here <f>ii and $/ 2 are the initial and final bound neu¬ 
tron states (for the definition of l\ and Z 2 see below). 
The cross section calculation result that the k 2 ~ ko = 
v / 2/ii 2 A/?i substitution may be used (see in Appendix 
D.) in calculating F\ and F 2 in a, where /ii 2 = too 

[(A\ — 1) (A 2 + 1 )] / (A\ + A 2 ). 

When evaluating (l23l) first the Weisskopf approxima¬ 
tion is applied, i.e. for the initial and final bound neu¬ 
tron states we take <&w (r„j) = (j>jw (r n j) Yi jmj (%), 
j = 1,2 where Yi jm . (Clj) is a spherical harmonics and 

<t>jw ( r nj) = yJs/Rj, J = 1,2 if |r nj -| < Rj and 

-i j o 

<t>jW (r n j) = 0 for \r n j\ > Rj, where Rj = roA- is 
the radius of a nucleus of nucleon number Aj with 
r 0 = 1.2 x 10 -13 cm. We apply the A x ~ A 2 ~ A\ — 1 ~ 
A 2 + 1 = A approximation further on. Calculating 
F\ (k 0 ) and F 2 (k 0 ) the long wavelength approximations 
(LWA) (exp (— *k 0 • r„i) = 1 and exp (?k 0 • r n2 ) = 1) are 
also used with s = 1 /ro that result approximately 


l 2 =rri2 

(|^(ko)| 2 ) Y, 1-^2 (ko)| 2 = 167r 2 ro/ 2 (2l 2 + 1). 
l 2 =—m,2 

(27) 

Using the results of Appendix D., the Ef i = A/2 
relation and if E e < 0.1 MeV (i.e. if F e (Ei e ) = 
\2irr] e {Ei e )\ = 2irZaf \Jm e (? /2£’j e ) then the cross sec¬ 
tion in the Weisskopf-LWA approximation reads as 


g w = — 


Cwo (2 1 2 + 1) r A2 A 3 ' 2 Z 2 


1 + 


2(A„—A) 


AA 


F e ( A/2) A 3 / 2 Et, 


(28) 


with Cwo = 2 9 7r 3 a 3 (0.08) 2 a^ro (^-) 3 (moc 2 ) 3 ^ 2 m e c 2 . 
Here as is the Bohr-radius, the relation c/v e = 
\Jm e c 2 / ( 2 Ei e ) with Ei e the kinetic energy of the in¬ 
going electrons is also applied and d = 3.52 x 10 “ 8 
cm (Ni lattice) and d = 3.89 x 10 -8 cm (Pd lattice). 
F e ( A/2) is determined by Q) and ©• The subscript 
W refers to the Weisskopf-LWA approximation and in 
(1251) the quantities A and E ie have to be substituted 
in MeV units. Cwo (Ni) = 1.4 x 10 -14 MeV 5 ^ 2 b and 
C W o (Pd) = 1.1 x lO " 14 MeV 5 / 2 b. 

We have calculated E; 2 =-m 2 \ F 2 (k 0 )| 2 , ^|Fi(k 0 )| 2 ^) 

and the cross section in the single particle shell model 
with isotropic harmonic oscillator potential and without 
the long wavelength approximation (see Appendix E.). 
We introduce the ratio 

(l^i(ko)l 2 ) g/i E(™ 2 |f 2 (ko)lL 

77 " ( i f >^) i 2 }^£-: 2 


(l^i(k 2 )| 2 ) 


1^2 (ko)l^ 


(29) 
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A 

58 

60 

61 

62 

64 

A_ 

-4.147 

-3.317 

0.251 

-2.526 

-1.587 

A+ 

0.928 

-0.251 

2.526 

-1.234 

-1.973 

r A 

0.68077 

0.26223 

0.0114 

0.03634 

0.00926 


TABLE I: Numerical data of the e + ^g 1 Ni + ^ 8 2 Ni —> e! + 
^ 1_1 IVi-|- ^ 2+1 Af» + A reaction. The reaction is energetically 
allowed if A = A_(Ai) + A+(A 2 ) > 0 holds. A is the mass 
number, ta is the relative natural abundance, A_ ( A ) = A a — 
Aa-i and A+(A) = Aa — Aa+i are given in MeV units. 


A 

102 

104 

105 

106 

108 

110 

A_ 

-2.497 

-1.912 

0.978 

-1.491 

-1.149 

-0.747 

A+ 

-0.446 

-0.978 

1.491 

-1.533 

-1.918 

-2.320 

ta 

0.0102 

0.1114 

0.2233 

0.2733 

0.2646 

0.1172 


TABLE II: Numerical data of the e + fg Pd +fg Pd —> e! + 
tg 1 Pd+fg +1 Pd + A reaction. The reaction is energetically 
allowed if A = A-(Ai) + A+(A 2 ) > 0 holds. A is the mass 
number, r A is the relative natural abundance, A_(A) = A a — 
Aa-i and A+(A) = Aa — Aa+i are given in MeV units. 


(The subscript Sh refers to the shell model.) With the 
aid of rj = ’nii,n 1 ,h,n 2 (Ai,A 2 ) given by (J93J) (see Ap¬ 
pendix E.) the cross section ash calculated in the shell 
model can be written as 

&Sh = Vh,ni,h,n 2 (Al,A 2 ) aw- (30) 


B. Yield of events of electron assisted neutron 
exchange process 


The yield dN/dt of events of electron assisted neutron 
exchange process A\,A 2 —> Ai — 1, A 2 + 1 can be written 
as 

dN 

— = N t N ni a$ , (31) 

dt 

where cr = {aw or ash }, $ is the flux of electrons, N t is 
the number of target particles, i.e. the number N Al of 
irradiated atoms of mass number Ai in the metal. The 
contribution of N n i neutrons in each nucleus X is also 
taken into account. N n i is the number of neutrons in 
the uppermost energy level of the initial nucleus ^fX. If 
F and D are the irradiated surface and the width of the 
sample, respectively, then the number of elementary cells 
N c in the sample is N c = FD/v c = 4 FD/d 3 in the case of 
Ni and Pd , and the number of atoms in the elementary 
cell is 2 rA x with r Al the relative natural abundance of 
atoms ^AX thus the number N f of target atoms of mass 
number A\ in the process is 

N t = J r Al FD . (32) 


The wave numbers and energies of the two outgoing 
heavy particles are approximately ki = —k 2 , 


E 1 


A 2 + 1 
A\ + A 2 


A and E 2 = 


Ai —1 

A\ + A 2 


A. 


(33) 


C. Numerical data of electron assisted neutron 
exchange processes in Ni and Pd 

As a first example we take Ni as target material. In 
this case the possible processes are 

e+ 2 £Ni + 28 Ni -+ e' + 2 i~ 1 Ni+ 2 g +1 Ni + A. (34) 


Ai ->• Ai - 1 

A .2 —y A 2 + 1 

A (MeV) 

V 

61 -> 60 

58 -> 59 

1.179 

7.02 x 10~ 3 

61 ->• 60 

61 -> 62 

2.777 

2.42 x 10" 8 

64—>63 

61 -> 62 

0.939 

2.08 x 10~ 4 


TABLE III: The values of the quantities r/ and A = A_ (Ai) + 
A-f (A 2 ) > 0, the later in MeV units, of the e + fNi + 
Ni — » e! + ^ 1_1 W-|- ^ 2+1 Ni + A reaction. TheA_(Ai) 
and A+(A 2 ) values can be found in Table I. For the definition 
of r] see (pll and 1(93)1. 


Tables I. and III. contain the relevant data for reaction 
m. Describing neutrons in the uppermost energy level 
of jg Ni isotopes we used lp shell model states in the cases 
of A = 58 — 60 and Of shell model states in the cases of 
A = 61 - 64. 

Another interesting target material is Pd in which the 
electron assisted neutron exchange processes are the 

e+ f£Pd+ 4 6 2 Pd —)• e' + fi^PdF f 6 2+1 Pd+ A (35) 

reactions. The relevant data can be found in Tables 
II. and IV.. Describing neutrons in the uppermost 
energy level of f e Pd isotopes we used 0 g shell model 
states in the cases of A = 102 — 104 and Id shell 
model states in the cases of A = 105 — 108. The nu¬ 
clear data to the Tables are taken from 0- One can 
see from Tables III. and IV. that in both cases three 
possible pairs of isotopes exist which are energetically 
allowed (for which A > 0) and their rates differ in 
the factor (2Z 2 + 1) N ni r} lltnit i 2tn2 (Ai, A 2 ) r Al r A2 A -3 / 2 
only. The rj = pq,n 1> j 2 , r t 2 (Ai, A 2 ) values of Ni and 


Ai —> Ai — 1 

A2 — y A2 ~h 1 

A (MeV) 

V 

105 -> 104 

102 -4 103 

0.532 

1.84 x 10 -4 

105 -> 104 

105 -> 106 

2.469 

8.88 x 10' 11 

108 -> 107 

105 —> 106 

0.342 

2.82 x 10 -3 


TABLE IV: The values of the quantities r/ and A = A_(Ai) + 
A + (A 2 ) > 0, the later in MeV units, of the e + ffPd +gg 
Pd — > e' + 4 g~ 1 Pd+fg +1 Pd+A reaction. The A_(Ai) and 
A+ (A 2 ) values can be found in Table II. For the definition of 
77 see dm and dS3]). 
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Pd can also be found in Tables III. and IV., re¬ 
spectively. The results of numerical investigation of 
(2 l 2 + 1 ) N ni rn uniMin2 (A 1 ,A 2 )r Al rA 2 A ~ 3/2 shows that 
the 61 —> 60, 58 —> 59 and the 108 —> 107,105 —> 106 
reactions are the dominant among the processes in Ni 
and Pd, respectively. 

In the case of Ni it is found that the 

e+ Ni+ jsfJVt -A e'+ HNi+ 5 2 9 s Ni+ 1.179 MeV (36) 

process of ash = 0.088 /E ie fib with E ie in MeV is lead¬ 
ing. In this case the 2 9 Ni and the 2 9 Ni isotopes take 
away 0.585 MeV and 0.594 MeV, respectively. In the 
case of Pd the 

e + \% 8 Pd+ \ 95 Pd -A e + ^Pd+ Pd + 0.342 MeV 

(37) 

reaction of ash = 0.26 /Ei e fib with Ei e in MeV is found 
to be the leading one. In this case the \ 97 Pd and the 
4 96 Pd isotopes take away 0.170 MeV and 0.172 MeV, 
respectively. 


III. OTHER RESULTS - OTHER CHARGED 
PARTICLE ASSISTED REACTIONS 

The transition probability per unit time and the cross 
section of the processes, which will be discussed below, 
may be determined in similar manner as was done above 
in the case of electron assisted neutron exchange process. 
The main difference is that in matrix elements V^ b and 
Vf* different Coulomb factors appear according to the 
particles which take part in the reaction. 


A. Electron assisted heavy charged particle 
exchange process 


There is an other possibility in the family of electron 
assisted exchange processes, when a charged heavy par¬ 
ticle (such as p, d, t, \He and 2 He) is exchanged. The 
process is called electron assisted heavy charged particle 
exchange process and it can be visualized with the aid of 
Fig.l too. Denoting the intermediate particle (particle 3 
in Fig. 1) by f 3 w, which is exchanged, the general elec¬ 
tron assisted heavy charged particle exchange processes 
reads as 


A i 


A' + 


z x -T z,Y^e' + %l£x*+ £+£y* + A. (38) 


Here e and e' denote electron and A is the energy of the 
reaction, i.e. the difference between the rest energies of 

initial {^\X +^ 2 v'j and final 
states. A = A. 


(fc 


3 X* 


A_i_ = Al 2 - A 


-^ 2+^.3 


A+, with A_ = A£ - A^i; 


+-23 
k A1—A3 


Y 


■) 


and 


, Ai 2 , Ai 2 +t are 


+ L ~*Z 2 L ~ 3 Z 2 +Z 2 ' L ~ X Z\ 5 ^Z\—Z 2 > t ~ X Z 2 1 L - i Z 2 -\-Z3 

the energy excesses of neutral atoms of mass number- 
charge number pairs A x , A x — A 3 , Z x — Z 3 ; A 2 , Z 2 , 
A 2 + A 3 , Z 2 + z 3 , respectively jlQ| . 



FIG. 2: The graph of electron assisted nuclear capture reac¬ 
tions. The simple lines represent free (initial (1) and final (1’)) 
electrons. The doubled lines represent free, heavy, charged 
initial (2) particles (such as p, d), their intermediate state 
(2’), target nuclei (3) and reaction product (4). The filled 
dot denotes Coulomb-interaction and the open circle denotes 
nuclear (strong) interaction. 


In (1551) the electron (particle 1) Coulomb interacts with 
the nucleus v)) X (particle 2). A scattered electron (par¬ 
ticle 1'), the intermediate particle (particle 3) and 
the nucleus 'z' l Zt^X* (particle 2 ') are created due to this 
interaction. The intermediate particle (particle 3) 
is captured due to the strong interaction by the nucleus 
f) 2 V (particle 4) forming the nucleus § 2 +^ 3 V* (particle 

5) in this manner. So in ( 1551 ) the nucleus ^X (parti¬ 
cle 2 ) looses a particle which is taken up by the 
nucleus f) 2 V (particle 4). The process is energetically 
forbidden if A < 0. Since particles 2', 3 and 4 all have 
positive charge, furthermore they all are heavy, the two 
Coulomb factors, which appear in the cross section, are 
F 2 ’ 3 and F34. Therefore the cross section of process ( 1551 ) 
is expected to be much smaller than the cross section of 
process ©• However process ( 1551 ) may play an essential 
role in explaining nuclear transmutations stated [I] (see 
below). Since Coulomb factors F 2 < 3 and F 34 determine 
the order of magnitude of the cross section of the pro¬ 
cess (the cross section of the process is proportional to 
F 2 > 3 F 33 ) we treat them in more detail in Appendix F. 


B. Electron assisted nuclear capture process 

Now the electron assisted nuclear caption process (see 
Fig. 2) is considered, in which an electron-nucleus 
Coulomb scattering is followed by a capture process gov¬ 
erned by strong interaction [g]. When describing the ef¬ 
fect of the Coulomb interaction between the nucleus of 
charge number Z and a slow electron one can also use the 
Coulomb factor F e (E) (0 of the electron defined above. 

As an example we consider the electron assisted d+d -A 
2 He process with slow deuterons. In this case, one of the 
slow deuterons (as particle 2) can enter into Coulomb 
interaction with a quasi-free, slow electron (as particle 
1) of the solid (see Fig. 2). In Coulomb scattering 
of free deuterons and electrons the wave number vec¬ 
tor (momentum) is preserved since their relative mo¬ 
tion may be described by a plane wave which is mul- 
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FIG. 3: The graph of heavy particle assisted nuclear capture 
reactions. The lines 1, 1’ represent free (initial (1) and fi¬ 
nal (T)) heavy particle which assists the reaction. The other 
lines represent heavy, charged initial (2) particles, their in¬ 
termediate state (2’), target nuclei (3) and reaction product 
(4). The filled dot denotes Coulomb-interaction and the open 
circle denotes nuclear (strong) interaction. 


tiplied by the corresponding Coulomb factor. In this sec¬ 
ond order process the Coulomb interaction is followed by 
strong interaction, which induces a nuclear capture pro¬ 
cess. The energy A of the nuclear reaction is divided be¬ 
tween the electron and the heavy nuclear product. Since 
toat m e (rriN is the rest mass of the nuclear product), 
the electron will take almost all the total nuclear reac¬ 
tion energy A away (there is no gamma emission) and 
the magnitude ky of its wave number vector k]/ reads 
ky = \/A 2 + 2m e c 2 A/ ( He ) ~ A/ (he) (if A ^$> m e c 2 ). If 
initially the electron and the deuteron move slowly and 
the magnitudes of their wave number vectors are much 
smaller than A/ (he), then the initial wave number vec¬ 
tors can be neglected in the wave number vector (momen¬ 
tum) conservation and consequently, in the intermediate 
state (in state 2') the deuteron gets a wave number vec¬ 
tor k 2 ' = — kp. If A = 23.84 MeV, which is the reaction 
energy of the d + d —> 2 F[e reaction, then the deuteron 2' 
will have k 2 ’ = A/ (he) and its corresponding (virtual) ki¬ 
netic energy Ey = A 2 / (4tooc 2 ) = 76.5 keV in the CM 
coordinate system. At this energy the Coulomb factor 
value between particles 2'and 3 reads as F 2 / 3 = 0.103. 
It must be compared to the extremely small Coulomb 
factor value, e.g. in the case of energy E = 1 eV to 
7*23 (1 eV) = 1.1 x 10 -427 , that is characteristic of the 
usual, first order process. If one compares again the cross 
sections of second order and first order (electron assisted 
and usual) processes then their ratio is approximately 
proportional to F e F 2 > 3 /F 23 (E) that becomes extremely 
large with decreasing E too. (The model and the de¬ 
tails of calculation, the results and their discussion can 
be found in Q.) The cross section of the electron assisted 
neutron exchange process is expected to be larger than 
the cross section of electron assisted nuclear capture pro¬ 
cess because of the appearance of the Coulomb factor in 
it. 


FIG. 4: The graph of heavy particle assisted heavy charged 
particle (such as p, d, t, \He and 2 He) exchange reaction. 
The lines 1, 1’ represent free (initial (1) and final (1’)) heavy 
particle which assists the reaction. The other lines repre¬ 
sent heavy, charged initial nuclei (2), their final state (2’, 
which is a nucleus lost particle 3), the transferred particle 
(3), target nuclei (4) and reaction product (5). The filled 
dot denotes Coulomb-interaction and the open circle denotes 
nuclear (strong) interaction. 


C. Heavy particle assisted nuclear processes 

In electron assisted nuclear reactions heavy, charged 
particles of energy of a few MeV may be created. In the 
decelerating process of reaction products of the electron 
assisted processes the energy of these heavy particles may 
become intermediately low (of about 0.01 [MeV]) so their 
Coulomb factor, if the particles are light, may be inter¬ 
mediately small so their assistance in nuclear processes 
have to be also considered among the accountable nuclear 
processes. The corresponding graphs can be seen in Fig. 

3 and Fig. 4. Fig. 3 depicts a heavy, charged particle 
assisted nuclear capture process and Fig. 4 represents 
heavy, charged particle assisted heavy charged particle 
(such as p , d, t, 2 He and 2 He) exchange reaction. Now 
all particles are heavy. According to the applied nota¬ 
tion, particles 2', 3 (in Fig. 3) and particles 3, 4 (in Fig. 
4) take part in a nuclear process and particle 1 only as¬ 
sists it. The different processes will be distinguished by 
the type of the assisting particle and also by the type of 
the nuclear process. In our model charged, heavy parti¬ 
cles, such as protons ( p ), deuterons (d) may be particle 
1 , which are supposed to move freely in a solid (e.g. in 
a metal). The other particles, that may take part in the 
processes are: localized heavy, charged particles (bound, 
localized p, d and other nuclei) as the participants of 
Coulomb scattering (with particle 1) and localized heavy, 
charged particles (bound, localized p, d and other nuclei) 
as nuclear targets (as particle 3 in Fig. 3 and particle 

4 in Fig. 4). The problem, that there may be identical 
particles in the system that are indistinguishable, is also 
disregarded here. 

The calculation of the transition probability per unit 
time of the process can be performed through similar 
steps to those applied for the calculation of the rate of 
an electron assisted process. The main difference is that 
now particle 1 is heavy. In order to show the capability of 
the heavy particle assisted nuclear processes, some cases 






of the proton assisted proton captures 

p+ + py 2 + 1 Y + p'+ A. ( 39 ) 

were investigated in Appendix III. (Cli. IX.) of ' 6 ). 

IV. DISCUSSION - ANALYSIS OF 
EXPERIMENTAL OBSERVATIONS 

A. Fleischmann-Pons experiment 

In the experiment of [l| Pd was filled with deuterons 
during electrolysis. The electrolyte had LiOD content 
too. Two types of electron assisted neutron exchange 
processes with Pd nuclei are possible: 

e + d + 4 ^ Pd —y e' + p Pd + A (40) 

with A = A —(d) + A+(A) and 

e +3 Li Pd —ye +3 Li + 4 ^"^ Pd + A (41) 

with A = A ~(Li) + A + (A) (the A+(A) values can be 
found in Table III). A_(d) = A^ — A p = 5.847 MeV 
and A -(Li) = A(jLi) — A(® Li) = 0.821 MeV are the 
energies of neutron loss of d and lLi, where A^, A p , 
A(lLi) and A(®Li) are the mass excesses of deuteron, 
proton, \Li and respectively. In reactions CU1) and 
m electrons of the metal are particle 1, d and 3 Li are 
particle 2 and f 6 Pd appears as particle 4 (see Fig. 1). 
Reaction (l40l) is energetically allowed for all the natural 
isotopes of Pd since A = A_(d) + A + (A) > 0 for each A 
(see the A + (A) values of Table II). In the case of reaction 
m the A = A-(Li) + A + (A) > 0 condition holds at 
A = 102 and A = 105 resulting A = 0.375 MeV and 
A = 2.312 MeV , respectively. 

However, at the Pd surface other types of electron as¬ 
sisted neutron exchange processes with d and Li nuclei 
of the electrolyte and d solved in Pd are possible: 


e + d d — y g! H - p t -|- A, 

(42) 

6 -b d -b d — y g! -b Tl -|- 2 He -b A, 

(43) 

6 ~b d H - \Li — y 6 -b p ~b 3 -^ ~b 

(44) 

e -b d ~b 3 Li — y c ~b 2 ^ He -b A, 

(45) 

6 -b d ~b 3 Li — y -b 2^ He -b Tl ~b A 

(46) 

e -b d -b 3 Li — y g! -b 4 Be -b 71 -b A, 

(47) 


which is promptly followed by the decay \Be —y 2\He 
(r Q = 6.8 eV). In these reactions electrons of the metal 
are particle 1 and d is particle 2 . 


In reaction (l40l) protons of energy up to 7.269 MeV 
and in reaction m 1 Li particles of maximum energy 
2.189 MeV are created which may enter into usual nu¬ 
clear reactions with the nuclei of deuteron loaded Pd and 
electrolyte which are (without completeness): the usual 
pd —y 2 -ffe + 7 reaction, 

p +3 Li -y 2%He + Q with Q = A + E kin (p), (48) 
®Li + d-y 2\He + Q with Q = A + E kin (Li), (49) 

3 Li T d — y p T 3 Li T Q with Q = A -t- E k i n (Li) . (50) 

In (1151) and (1151) the emitted fife has energy Ei He > 
8.674 MeV and E* He > 11.186 MeV, and in (l50l) the 
created p and 3 Li have energy E p > 4.397 MeV and 
Er Li > 0.628 MeV, respectively. It can be seen that in 
(1151) and (1151) 2 -ffe is produced. The 3 Li particles may 
enter into reaction 

3 L 2 T d —^ 2 2 ffe T Ti T Q with Q = A -f- E k i n (Li) (51) 

which contributes to the \He production too. Here and 
above E k i n (p) and E k i n (Li) are the kinetic energies of 
the initial protons, 3 Li and 3 Li isotopes. 

From the above one can see that at least twelve types 
of reactions (altogether 18 reactions) exist which are ca¬ 
pable of energy production and in half of them energy 
production is accompanied with \He production. It is 
reasonable that reactions (HU1) and (1H1) have the high¬ 
est rate in the above list of reactions. In the majority 
of the above reactions charged particles, mostly heavy 
charged particles are created with short range and so 
they loose their energy in the matter of the experimental 
apparatus mainly in the electrode (cathode) and the elec¬ 
trolyte, therefore their direct observation is difficult. It is 
mainly heat production, which is a consequence of decel¬ 
eration in the matter of the apparatus, that can be expe¬ 
rienced. The third of the processes, mainly the secondary 
processes are the sources of neutron emission. X- and 
7 —rays may be originated mainly from bremsstrahlung. 
The above reasoning tallies with experimental observa¬ 
tions. 

In reactions (|551) — (l5ll) heavy, charged particles of ki¬ 
netic energy lying in the MeV range are created which 
are able to assist nuclear reactions. One can obtain the 
possible heavy charged particles assisted reactions if in 
reactions m - m the electron is replaced by heavy 
charged particles (p, t, |ffe, \Ee, %Li, \Li, \Be and 
^ 6 +1 Pd with A = 102,104 — 106,108,110) which are cre¬ 
ated in reactions flUl) — (1511) . Since the number of possible 
heavy charged particles is 13 and the number of reac¬ 
tions which may be assisted by them is 8, at least 104 
heavy charged particle assisted reactions must be taken 
into account. Consequently, it is a rather great theoreti¬ 
cal challenge and task to determine precisely the relative 
rates and their couplings of all the accountable reactions, 
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a work which is, nevertheless, necessary for the accurate 
quantitative analysis of experiments. 

The relative rates of coupled reactions of many types 
depend significantly on the geometry, the kind of matter 
and other parameters of the experimental apparatus and 
on some further variables, which may be attached to a 
concrete experiment. This situation may be responsible 
for the diversity of the results of experiments, which are 
thought to have been carried out with seemingly in the 
same circumstances. 


B. Nuclear transmutation 

As to the phenomenon of nuclear transmutation [|j we 
demonstrate its possibility only. First let us see the possi¬ 
bility of normal reactions. For instance in a Fleischmann- 
type experiment ® Li particles of energy up to 2.189 MeV 
are created in reaction dSD so the reaction 

3 Li +3 Li — C + 'y + Q (52) 

may have minor, but measurable probability. Here Q = 
A + Ekin(Li). 

The Coulomb factor of reaction (15^1) is FLi,Li = 1-71 x 
10 -3 at 2.189 MeV kinetic energy of %Li particles. The 
magnitude of the Coulomb factor indicates that the rate 
of reaction (15^1) may be large enough to be able to pro¬ 
duce carbon traces in observable quantity. 

Moreover, in reactions (HU1) and (PHI) free f e Pd particles 
are created offering e.g. the possibility of the 

e Pd Pd — y e! +44 4 Ru +43 Cd - 1 - A (53) 

electron assisted ^He exchange process. The electron and 
the other Pd particle are in the solid. Analyzing mass 
excess data [lC| it was found that e.g. the e+\g 3 Pd+\^ 
Pd —y e' +44° Ru +4I 4 Cd + A \FLe exchange process has 
reaction energy A = 5.7305 MeV. \^Pd and l^Pd are 
produced in reaction (l40l) .] Calculating the TV 3 = F34 
Coulomb factors taking A = 100, Z = 46, A 3 = 3, Z 3 = 2 
in (1103D one gets FVs-Fm = 2.5 x 10 -12 which seems to 
be large enough number to produce Cd and Ru traces in 
an experiment lasting many days long. 

The above reactions may offer starting point for the 
explanation of nuclear transmutations. 

C. Rossi-type reactor (E-cat) 

Recently the Rossi-type reactor [111 (E-Cat) was ex¬ 
perimentally investigated in detail 12]. The fuel con¬ 
tained mostly Ni and also Li in accountable measure, 
there was O.Ollg Li in 1 g fuel. The isotope composi¬ 
tion of the unused fuel was equal to the relative natural 
abundances. But the isotope composition of the ash (the 
fuel after 32 day run of the reactor) strongly changed. 
(The measured relative abundances of Li and Ni iso¬ 
topes in fuel and ash can be seen in Table V. The natu¬ 
ral abundances are also given for comparison. The data 


I sotope Fuel 

Ash 

Natural 

1 Li 

0.086 

0.921 

0.075 

3 Li 

0.914 

0.079 

0.925 

i m 

0.67 

0.008 

0.681 

1 ° s Ni 

0.263 

0.005 

0.262 

as Ni 

0.019 

0.000 

0.018 

28 Ni 

0.039 

0.987 

0.036 

28 Ni 

0.01 

0 

0.009 


TABLE V: Measured relative abundances of Li and Ni iso¬ 
topes in fuel and ash. The natural relative abundances are 
also given for comparison. The data are taken from fl2[. 

are taken from Appendix 3. of 0.) One can see that 
the 2 g./V 7 isotope is enriched and the other Ni isotopes 
are depleted. Furthermore, the relative 3 Li content de¬ 
creased from 0.917 to 0.079 while the relative %Li content 
increased from 0.086 to 0.921. 

The reactor worked for about ten days at temperature 
T\ = 1533 K and the remaining time at temperature 
X 2 = 1673 K. At these temperatures a free electron gas 
may be created from the Ni powder of the fuel due to 
the termionic emission process. The emitted flux of elec¬ 
trons can be determined from the current density of elec¬ 
trons according to the Richardson’s law using the work 
function U = 5.24 eV of Ni. The obtained termionic 
electron fluxes are $1 = 7.5 x 10 9 cm~ 2 s~ 1 and $2 = 
2.4 x 10 11 cm~ 2 s~ 1 at Ti and T 2 , respectively. Regard¬ 
ing the large surface of the powder fuel it is reasonable 
to suppose that the free electron gas is formed near the 
surfaces of grains of powder. But if a free electron gas 
interacts with the LiAlH± — Ni powder mixture applied 
then the above observations can be well explained by the 
electron assisted neutron exchange processes. 3 Li has 
A_ = 0.8214 MeV so it is able to lose neutron. The 
A + values of the Ni isotopes can be found in Table I. 
Completing Table I with A_|_(!|./V*) = 3.319 MeV (the 
half life of | %Ni is r = 7.6 x 10 4 y) one can recognize that 
the e+ \Li- V —y e'+ 3 Li+ ^ 8 + 1 A r * + A reaction has 

A > 0 value for A = 58 — 61 but in the case of A = 62 the 
chain of reactions breaks since in this case A < 0 because 
A+(%Ni) = -1.234 MeV. The 64 63; 61 -> 62 reac¬ 

tion of type (1M1) (see Table III) leads to production of 
2 glV* (r = 100.1 y) which has A-dglV*) = 1.2335 MeV 
allowing and coupling transition 63 —y 62 to transitions 
58 —y 59; 59 —y 60; 60 —y 61 and 61 —> 62 in reaction (f34|) . 
These facts explain the enrichment of IgA’* and 3 Li and 
the depletion of 3 Li and Ni isotopes of A = 58 — 61 and 
64. (Reactions (l34l) too contribute to the enrichment of 
2 \Ni (see Table III).) 


V. CONCLUSION 

It is thought that, in principle, the electron assisted 
processes are able to answer the questions raised in the 
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introduction. The exchange of the original, extremely 
small Coulomb factor to the Coulomb factor of order of 
unity of the electron in electron assisted processes an¬ 
swers problem (a). The electron assisted nuclear reac¬ 
tions and the reactions which are coupled with them are 
not accompanied by the expected nuclear end products 
answering problem (b). Problem (c), the asserted ap¬ 
pearance of nuclear transmutations is partly answered in 
Section IV.C. with the aid of charged particle assisted 
and usual nuclear reactions. 

Summarizing, the theoretical results expounded and 
their successful applications in explaining some unre¬ 
solved experimental facts inspire us to say that the study¬ 
ing of charged particles electron assisted nuclear reac¬ 
tions, especially the electron assisted neutron exchange 
processes may start a renaissance in the field of low en¬ 
ergy nuclear physics. 


VI. APPENDIX 

A. Initial, intermediate and final states of electron 
assisted neutron exchange process 

Let H/j, and 'L/ denote the space dependent parts 
of initial, intermediate and final states, respectively. The 
initial state has the form 

$ ! (x e ,Xi,X n i,X 2 ) = i> ie (x e ) 1 p iln (xi, X„i) 1 p i2 (x 2 ), 

(54) 

where 


1pie (X e ) = p-l/2 e fik ie -x e ) and ^. 2 ( X2 ) = y-l/2 e (ik i2 .x 2 ) 

(55) 

are the initial state of the electron and the nucleus ^ 2 X, 
and ^iin(xi, x„i) is the initial state of the neutron and 
the initial A\ — 1 nucleon of the nucleus ^X. x e , x x , x nX 
and x 2 are the coordinates of the electron, the center of 
mass of the initial A\ — 1 nucleon, the neutron and the nu¬ 
cleus f) 2 X 1 respectively, k i e and k; 2 are the initial wave 
vectors of the electron and the nucleus ^ 2 X and V is the 
volume of normalization. The initial state V , *in(xi, x„ x ) 
of the neutron and the initial A\ — 1 nucleon may be given 
in the variables R x , r„ x 

Viii„(Ri, r nX ) = V -1 / 2 exp(*k,i • R x )3> ix (r„ x ) (56) 

where Ri is the center of mass coordinate of the nu¬ 
cleus z*X and r„ x is the relative coordinate of one of 
its neutrons. R x and r„ x are determined by the usual 
x„i = Ri + r„i and Ri = [(A x - 1) x x + x„ x ] /A ± rela¬ 
tions where x„ x and x x are the coordinates of the neu¬ 
tron and of the center of mass of the initial A 1 — 1 
nucleon, respectively. The inverse formula for x x is 
xi = Ri — r„i/ (Ai - 1). In (HU) the 4> ix (r„ x ) is the 
wave function of the neutron in the initial bound state 
of nucleus ^CY, k^is the initial wave vector of nucleus 

z 1 *- 


The intermediate state has the form 
4' M (x e ,Xi,X„i,X 2 ) = Xpf e (x e )'i/> At i n (xi,X n i)V’i2(x 2 ), 

(57) 

where 

Ipfe (X e ) — 1 / 2 gb k /e -x e) (5 8 ) 

with k f e the wave vector of the electron in the final state 
and ipi 2 (x 2 ) is given in (l55l) . The state x„ x ) is 

the product of two plane waves ipf x (x x ) = V _ 1 / 2 e( zk 1 ' Xl ^ 
and ip n (x„i) = V'- 1 / 2 e * k »- x »i ) which are the final state of 
the nucleus J ^~ 1 X and the state of the free, intermediate 
neutron. Thus Vvin(xi, x„ x ) = V - 1 e 8k i x i e lk «' x "i and it 
has the form in the coordinates R x , r nX 



where k x and k„ are the wave vectors of the nucleus 
and the neutron, respectively. 

The intermediate state may have an other form 


’JV(x e ,Xi,X„i,X 2 ) = l/lfe (x e )^/i(Xi)VV2n(Xnl,X 2 ), 

(60) 

where 


1 p^ 2 n(,X n 1 ,X 2 ) = lj) n (x n l )^* 2 (x 2 ) = V V krl ' Xnl e * k i 2 -x 2 

( 61 ) 

which can be written in the coordinates r n2 = x„ x — R 2 
and R 2 = (A 2 x 2 + x nl ) / (A 2 + 1) as 


n (R 2 7 Y n 2 ) 


J_„d k i 2 + k „)'R 2 „ l ( k " 

V 



(62) 


where R 2 is the center of mass coordinate of the nucleus 
z 2+1 X and r„ 2 is the relative coordinate of the neutron 
in it. In these new variables x 2 = R 2 — r n 2 /A 2 and 
x„ x —x 2 = (A 2 + 1) r n2 / A 2 which is used in the argument 
of V St (given by (fl2l) l in calculating Vj*. Evaluating the 
matrix elements V^ b and Vj* the forms (15U1) and (linil) of 
tpfj, are used, respectively, and d 3 k n in ED- 

The final state has the form 


^/(Xe, Xi, X„i, X 2 ) = Ip fe (x e ) Ipfl (xi)lpf2n (x nX , X 2 ), 

(63) 

where ^/ 2n (x„i,x 2 ) is given in the variables R 2 , r n2 as 

V’/ 2 n(R 2 , r n2 ) = V ~ 1/2 exp(ik 2 • R 2 )$/ 2 (r n2 ) , (64) 

and $/ 2 (r n2 ) is the bound state of the neutron in the 
nucleus ^ 2 + 1 -Y. 


B. Evaluation of matrix elements V^ b and VfP 


The argument of the Coulomb potential V Cb is x e — x x 
therefore the integration with respect to the components 
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of x 2 may be carried out and / |'0z2 (^ 2 ) | 2 d 3 £ 2 = 1. The 
remainder is 

V™ = J V>/ e ( x e)^l„(xi,X n i)y C ' 6 (x e -Xi) (65) 

xip ie (x e ) ipii n (xi,x. n i)d 3 x e d 3 xid 3 x n i. 

Making the X]. x n] —> R 1; r„i change in the variables, 
substituting the forms (l56l) and (l59l) of ifn n and ip^xn, 
and neglecting k,i, the integrations over the compo¬ 
nents of x e and Ri result V~ 3 (27r) 3 5 (q + k ie — kf e ) 
and V~ 3 / 2 {2tt) 3 S (q — ki — k n ), respectively and the in¬ 
tegration over the components of r nl produces Fx (k„) 
where 

F\ (k„) = j ( h,;i (r„i) e _ *( k,1_A i- 1 )' r ” 1 d 3 r nl . (66) 


Using the 6 (q + k ie — k/ e ) in carrying out the integra¬ 
tion over the components of q in V^ b one gets 


V Cb 

r /ll 


Ze 2 


27r 2 |k fe - k ie | 
x \/G^5 (k ie ^ k fe - ki - k 


-F, (k„) <0 x (67) 


and 


Fx (k„) 


_ / k l+ k /e- k ie \ 

i(r„i)e ^ > nl d 3 r nl . 


( 68 ) 

For particles e and 1 (ingoing electron of charge —e and 
initial nucleus ^fX of charge Ze) taking part in Coulomb 
interaction we have used plane waves therefore the matrix 
element must be corrected with the so called Sommerfeld 
factor 0 a/Gs where 


F e (Ei e ) 
~ Fe(Ef l)' 


(69) 


Now we deal with Vf*. The strong interaction works 
between the neutron and the nucleons of the nucleus X 
therefore the argument of V St is x„i — x 2 . The integra¬ 
tions with respect to the components of x e and xi result 
f \ipef(*-e )\ 2 d 3 x e = f |^/i(xi)| 2 d 3 Xi = 1. The remain¬ 
der is 



x-2)i’,i2nd 3 X2d 3 x nl . 


(70) 


Similarly to the above, making the x„i,x 2 —> R 2 ,r n2 
change in the variables, substituting the forms (l62l) and 
(1M1) of and ^f 2 n an d neglecting k i2 , the integrations 

over the components of R 2 result V~ 3 / 2 (27t) 3 S (k„ — k 2 ) 
and the integrations with respect to the components of 
v n 2 produces F 2 (k„) with 


F 2 ( k„) = 


$ 


/2 ( r « 2 )e 


ik„-r„ 2 


(71) 


-/ 


exp(-s^±ir„ 2 ' 

W ^2 , 


d T n 25 


where r n2 = |r„ 2 |. Taking into account that the neutron 
interacts with each nucleon of the final nucleus of nucleon 
number A 2 

vf* = 5 ( k « - k 2 ) • ( 72 ) 


C. Transition probability per unit time of electron 
assisted neutron exchange process 


Substituting the obtained forms of V^ b and Vf* (for¬ 
mulae (l67l) and (T72l) 1 into m and using the correspon¬ 
dence V -»■ jA n d 3 k n and the S (k„ — k 2 ) in the inte- 

z —'A 4 v ' 

gration over the components of k.„ one gets 


Tfi = - 


e 2 ZA 2 F 1 (k 2 )U 2 (k 2 )^fg^ 


27T 2 |k f e - k 


A 2 


(73) 


(27t) 6 (5 (ki + k 2 + kfe - k ^) 


U 3 


(A-E / i*) k _ k 


where 
Fx (k 2 ) = 


$ii (rnl) e 


■ ( , kl+k/e-kieN 

2 \ 2 - 4-1 ) 


r „l ,3 


drr n x (74) 


and F 2 (k 2 ) is determined by (l26l) . Here and 
<f> f 2 in (1261) are the initial and final bound neu¬ 
tron states. Substituting the above into m, 
using the identities [5 (ki + k 2 + kf e — k; e )] 2 = 
S (ki + k 2 + kfe — k^) d (0) and (27r) 3 <5 (0) = V, 


X 3 


the *■ Em I d 3 kxd 3 lt 2 d 3 kf e correspon¬ 


dence, averaging over the quantum number m 1 and 
integrating over the components of kf e (which gives 
kf e = —ki — k 2 + ki e ) one obtains 


W fi = 


1 


a}hc 2 Z 2 Y!C- 2 m 2 l^(k 2 )| 2 (?5) 

tt 3 v c V (|ki + k 2 | 2 + A 2 )" (A^)k„= k2 

( |Fl (kz)|2 ) |rjf^W(£/ ~ 


where Ax, A 2 are the initial atomic masses, h,mx and 
Z 2 ,to 2 are the orbit and its projection quantum numbers 
of the neutron in its initial and final state. For F\ (k 2 ), 

(\Fx (k 2 )| 2 ^) and F 2 (k 2 ) see (O, d25]) and ([251) . Taking 
into account the effect of the number of atoms of atomic 
number H 2 in the solid target the calculation is similar 
to the calculation of e.g. the coherent neutron scattering 
0 and the |7/i| 2 must be multiplied by Nl which is the 
number of atomic sites in the crystal and by r^ 2 which is 
the relative natural abundance of atoms ^ZX. We have 
used Nl/V = 2 /v c with v c the volume of the elementary 
cell of the fee lattice in which there are two lattice sites 
in the cases of Ni and Pd investigated. 
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D. Approximations, identities and relations in 
calculation of cross section 


E. (|Ei (k 0 )| 2 ) sh and £^=-^ 2 \F 2 (k 0 )|| h in single 
particle shell-model and without LWA 


Now we deal with the energy denominator (A E^i) in 
(1751) and (1751) [see (1l5l) — (1771) ]. The shielding parameter 
A is determined by the innermost electronic shell of the 
atom z x X and it can be determined as 


a = A 

as 


(76) 


where as = 0.53 x 10 8 cm is the Bohr-radius. The 
integrals in (1751) and (1751) have accountable contributions 
if 


|ki + k21 A 


(77) 


and then Ef e < h 2 A 2 / (2 m e ) = | a 2 m e c 2 Z 2 which can be 
neglected in A E^i and in the energy Dirac-delta. Thus 


A E^i — 


and in the Dirac-delta 


E t = 


! h 2 ko 

(78) 

■ + „ A_ + An 

ZiTYl n 

h 2 k? h 2 k 2 

(79) 

1 + 2 . 

2 mi 2 m 2 


In this case ki = — k 2 + (5k with |<5k| = Sk ~ A. Using 

ki ~ fc 2 ~ k 0 = \f 2 [i\ 2 A/h (80) 

(see below) with /ii 2 c 2 = Ai 2 tooc 2 , where Ai 2 = 
(Ai — 1 ) (A 2 + 1 ) / (Ai + A 2 ) is the reduced nucleon 
number, one can conclude that the k 2 = —ki relation 
fails with a very small error in the cases of events which 
fulfill condition (1771) since fci/fco — 1 , k 2 /ko — 1 , 5k/ko ~ 
X/k 0 and X/ko = a.f Zm e c 2 /\j2gi 2 c 2 A -C 1. Conse¬ 
quently, the quantity Ef in the argument of the energy 
Dirac-delta can be written approximately as 

Ff = (JL + JL\ k 2 = -^^_ fail 

f \ 2 mi + 2 m 2 ) 2 2 A 12 m 0 c 2 ' ^ ' 

Furthermore taking A\/ (Ai + 1) ~ 1 


h 2 c 2 k| 

— — -— — A_ + A n . 


2 moc 2 


(82) 


We introduce the Q = fick 2 /A, P = He ((5k) /A, £f = 
Ef /A = [Q 2 / (2Ai 2 moc 2 )] A and L = hcX/A dimen¬ 
sionless quantities. The energy Dirac-delta modifies as 
6 (Ef — A) = S[ef (Q) — 1] /A. The relation (1751) yields 
L = hcZ/ (as A) = Zajm e c 2 /A and Zctfm e c 2 /A < 1. 
Now we change d 3 k\d 3 k 2 to (J ^) 6 d 3 Qd 3 P in the inte¬ 
gration in ([75]), use the 5 [g (Q)] = 6 (Q — Q 0 ) /g'(Q 0 ) 
identity, where Q o is the root of the equation g ( Q ) = 0 
(ko = QqA/ (he), see (1501) 1. estimate the integral with 
respect to the components of P by 

f 00 4np2dP = (A (83) 

Jo (P 2 + L 2 ) 2 L 

and apply v c = d 3 / 4 (the volume of unit cell of fee lattice 
for Ni and Pd of lattice parameter d). 


Now we calculate the quantities ^|Ei(k 0 )| 2 ^ and 

Xu 2 =-m 2 1-^2 (ko)l 5 /, in the single particle shell model 
with isotropic harmonic oscillator potential and without 
the long wavelength approximation (see definitions: (Hfl), 
(1751) and ©). Taking into account the spin-orbit cou¬ 
pling in the level scheme the emerging neutron states are 
01 and 11 shell model states in the cases of Ni and Pd 
to be discussed numerically ITU . So the initial and final 
neutron states ($u, $/ 2 ) have the form 


*Sh (r nj ) = (ft.) 

‘nj 


(84) 


where rij = 0,1 in the cases of 0 Z and 1 / investigated, 
respectively, and 


Ro h = b ~ 1/2 ( 


2 h+l ( 1 2 

T(lj +3/2)) * 6XP \ 2^ J 


1/2 


(85) 


A'x (86) 
x 11 - WT 3 e /> exp (T? 


IS 


with Qj = r n j/bj where bj = y 'h/ ( moujj) |l5[. Here Wj 
the angular frequency of the oscillator that is determined 

by huii _= 40 A ^ 1 ^ 3 MeV and hui 2 = 40 (A 2 + l ) -1 ^ 3 

MeV [16|. (The subscript Sh refers to the shell model.) 
With the aid of these wave functions and for ni = 0,1 

2 Zl+2 


with 


< 87 > 

/*°° ^ 2 
h,o= Q h+2 jh( k obi . 1 g)e~i e2 dg (88) 
Jo A i “ 1 


and 


h,i = l 


Ji+2 


1 - 


2Zi + 3 C 


g 2 x (89) 


A i 


xjh(k 0 bi- - -g)e * e dg. 

A± — L 

Here ji x (x) = yf^Ji 1 + i/ 2 (x) denotes spherical Bessel 
function with Jq+ 1 / 2 (x) the Bessel function of first kind. 
Similarly 

= (90) 

l 2 ——m2 


X47T 


A 2 


a 2 + 1 


'2,712 
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with 


^ 2 , 0 = / Q l2+1 ji 2 {kob2Q)e 28 ~^~ r o e dg (91) 


and 


I2,1 — l 


0 h + 1 ( 1 - 


2l 2 + 3‘ 


p 1 2 3 4 x (92) 


xj; 2 (M 20 )e 28 r o e dg. 


Substituting the results of (l87l) . (l90l) and (l27l) into 
one gets 


Vh,1^1,12,112 (Al, A?) — 


2 ^ 1 +^ 2+4 


7r [2l\ + 1)!! (2l 2 + 1)!! 

b\b 2 f A 2 k 


x (93) 


A 2 + 1 


t2 t2 
1 l,n 1 1 2,n 2 ■ 


F. Coulomb factors F 2 > 3 and F 34 in electron assisted 
heavy charged particle exchange process 


Calculating the Coulomb factor F 2 i 3 [see ®] between 
particles 2' and 3 the energy determined by (l97l) is given 
in their CM coordinate system (since k 3 = — k 2 ') thus it 
can be substituted directly in fU producing 


V 2'3 = (Z 2 - £ 3 ) z 3 a f A 3 i 


A 2 ' + A 5 nioc 2 


(98) 


(A 2 / + A 3 ) A 5 2A 
Since the above analysis is made in order to discuss the 
phenomenon of nuclear transmutation we take A 3 -C 
A 2 i ~ A 3 = A (> 100 in the case of Pd discussed). 
So (A 2 i + A 5 ) / [(A 2 / + A 3 ) A 5 \ ~2/A and 77^3 reads ap¬ 
proximately as 


V2'3 — {Z 2 


z 3 ) z 3 otfA 3 



(99) 


Calculating the Coulomb factor F 34 , the energy of parti¬ 
cle 3 determined by (1971) is now given in the laboratory 
frame of reference since particle 4 is at rest. In the CM 
system of particles 3 and 4 the energy E 3 {CM) is 


If initial particles have negligible initial momentum 
then, because of momentum conservation, k 2 / = ks in 
the final state. (It was obtained 0 that the process has 
accountable cross section if the momentum of the final 
electron can be neglected, i.e. in the kp ~ 0 case.) Thus 
the condition of energy conservation 


^ 2 kn, fi 2 k? 

- — H- 2 . = a 

2 to 2 ' 2 to 5 


(94) 


determines k 2 / as 


h 2 k 2 , = 2g, 2 i 5 A , (95) 


where H is the reduced Planck-constant, 

M 2 ' 5 , = a 2 ' 5 ?TloC 2 (96) 


is the reduced rest mass of particles 2' and 5 of mass 
numbers A 2 / and A 5 [for a 2 ’ 3 see (|5|)]. If the initial 
momenta and the momentum of particle V are negligi¬ 
ble then k 3 = — k 2 /, since momentum is preserved in 
Coulomb scattering. Thus the energy E 3 of particle 3 
can be written as 


E 3 (CM) 


A 4 a 2 > 3 A 
(A 3 + A 4 ) A 3 


( 100 ) 


Substituting it into fH) 


7734 = (Z 4 + 2 : 3 ) z 3 a f A 3 


1 moc z 
2a 2 '5 A 


( 101 ) 


Applying the same approximation as above in which 
2a 2 '5 ~ A 


V 34 : — {Z 4 + z 3 ) z 3 afA 3 



( 102 ) 


Furthermore, if Z 2 ~ Z 4 = Z z 3 then 


r/ 2'3 = 7734 = Zz 3 a,fA 3 



(103) 


E 3 


h 2 H = P2X A = a 2'5 ^ 

2 m 3 m 3 A 3 


(97) consequently, F 2 , 3 = F 34 . 
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